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Abstract. We present mass distributions obtained from 
near-infrared (NIR) surface brightness decompositions 
and rotation curve fitting of a sample of early-type spi- 
ral galaxies. Bulge and disk mass-to-light (M/L) ratios, 
dark halo parameters, and the modified Newtonian dy- 
namics (MOND) acceleration parameter are derived. We 
find that the mean disk NIR M/L is higher than that of 
the bulge, and comparison with stellar population syn- 
thesis models implies that early-type spiral bulges are, on 
average, younger and more metal rich than disks. NIR 
disk M/L is found to depend on disk luminosity, consis- 
tently with previously reported trends in ellipticals and 
spirals, and with cold dark matter models for disk forma- 
tion. Within the optical radius, the mean ratio of stellar to 
dark matter is 2 and the typical dark halo mass is 10^^ Mq. 
The value of the MOND acceleration parameter that best 
accommodates the sample as a whole is 1.3 • 10~® cms"^. 

Key words: Galaxies: spirals - Galaxies: fundamental 
parameters - Galaxies: structure - Galaxies: kinematics 
and dynamics - Galaxies: stellar content 



1. Introduction 

The comparison of the luminosity profiles and rotation 
curves of disk galaxies provided the first evidence for dark 
mat ter in outer parts of galaxies (e.g., van Albada & San- 
1986 ). The luminous material, either in the form of 



cisi 



stars or neutral gas, is not able to reproduce the approxi- 
mate flatness of the rotation curves at large radii, implying 



that the Newtonian theory of gravitation needs to be mod- 
ified or that the mass in the outer regions of most spirals 
is dominated by a dark component. The physical extent 
of such a "dark halo" cannot be inferred from the rotation 
curves, since these last are flat even where the luminous 
material ceases to be detected. Hence, the total mass in 
spiral galaxies is not yet well known. 

The properties of dark matter (DM) and the structure 
of dark halos as derived from rotation curves have pro- 
found implications on cosmological issues such as galaxy 
formation. Recent N-body simulations of cold dark mat- 
ter halos together with adiabatic infall models have shown 
that matching observed rotation curves requires a system- 
atic increase of disk mass-to-light ratios with luminosity 
(Navarro et al. |1996D . Such a trend is observed in large 
bodies of rotation curve data (Broeils 1992| ; Persic et al. 
1996 - hereafter PSS), and further comparisons between 
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theoretical predictions and observations may help clarify 
the nature of dark matter in galaxies. 

The presence of a dark halo makes global mass-to-light 
(M/L) ratios derived from rotation curves much higher 
than those expected from stellar populations, and pro- 
duces radial gradients in the total M/L. The contribution 
of DM tends to be larger for low-luminosity systems, as 
the fraction of DM to luminous matter within the optical 
radius appears to decrease with incre asing galaxy lumi- 
nosity (Kormendy 199C; Salucci et al. 1991). In luminous 
systems, therefore, the mass distribution in the inner few 
disk scale lengths can be largely ascribed to the bulge 
and disk stars. Consequently, when stellar galaxy compo- 
nents can be accurately isolated with surface photome- 
try, estimates of their M/L from rotation curves provide 
useful constraints for the age, abundance, and star forma- 
tion history of their stellar populations. Recent population 
synthesis models predict trends of M/L with metallicity 
and with wavelength (Worthey 1994; Bruzual & Chariot 
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1993), and can be used to infer and compare properties of 
the stellar populations in bulges and disks. 

We have reported in a previous paper (Moriondo et al. 
hereafter Paper I) the results of two-dimensional 



1997 



near-infrared (NIR) surface brightness decompositions for 
a sample of early-type spirals. NIR wavelengths, especially 
when combined with those in the optical, are a powerful 
diagnostic tool since they trace more accurately the stellar 
mass content and minimize the complications of extinc- 
tion. In this paper, we apply the bulge-|-disk decomposi- 
tion results to the analysis of the rotation curves for our 
sample. In Sect. 2 we derive the radial mass distributions 
and evaluate the contribution of the bulge, disk, and dark 
halo to the observed rotation curves; as an alternative to 
dark matter halos, we also derive a value of the modi- 
fied Newtonian dynamics critical acceleration parameter. 
The resulting NIR M/L ratios of the luminous compo- 
nents are analyzed in Sect. 3, and compared with those at 
optical wavelengths and with population synthesis mod- 
els. Finally, we assess the importance of dark halos, and 
explore trends of M/L with luminosity in the context of 
the properties of the fundamental pl ane fo r gravitationally 
bound objects (e.g., Burstein et al. 1997 ). 



2. Data analysis 

2.1. Surface photometry and .structural decomposition 

The techniques used to reduce and analyze the photomet- 
ric data are detailed in Paper I. Briefly, we decomposed the 
surface brightness distributions of each sample galaxy into 
the two components of bulge and disk. Two independent 
methods were used: a) a parametric fit assuming a gen- 
eralized exponential bulge plus an exponential thin disk; 
and b) an iterative non-parametric (np) decomposition al- 
gorithm. We emphasize that both procedures decompose 
the entire two-dimensional brightness distribution, rather 
than the brightness profiles, and take into account the 
effects of seeing. The bulge is assumed to be an oblate 
rotation ellipsoid coaxial with the disk. In addition to the 
NIR images, r band brightness profiles and M/L ratios 
(Kent 1988) are available for ten galaxies of our sample. 



2.2. Observed rotation curves 

Rotation curves (RC's) exist for all the galaxies in our 
sample; the references for them are shown in Table 0, 
together with general information regarding the sample. 
More information about the sample objects (coordinates, 
etc.) can be found in Paper I. HI data are available for six 
objects; for the remainder, RC's were obtained from opti- 
cal emission lines. We discarded absorption data since they 
are likely to yield a lower limit rather than a true estimate 
of the rotational velocity, due to the effect of velocity dis- 
persion and projection along the line of sight. All distances 



were scaled to match our values, and circular speeds were 
corrected when estimated from inclinations differing from 
ours, with the exception of extended HI RC's derived from 
interferometric maps. 

2.3. Model rotation curves 

To determine the shape of the RC for our galaxies we as- 
sumed: a) optical transparency (which was also assumed 
for the surface brightness decompositions); and b) the 
M/L ratio of bulge and disk to be constant within each 
component. The circular velocity in the equatorial plane 
of an axisymmetric mass distribution is given by 



w^(r, 0) = r 



9r ) (z=o) 



(1) 



where $(r) is the gravitational energy at radius r. In this 
case $ is the sum of the contributions of bulge and disk 
and, when necessary, also of a dark halo. 

2.3.1. The bulge 

The circular speed of an axisimmetric ellipsoida l dist ribu- 
tion of matter is given by (Binney & Tremaine 1987): 



vl{r)=4TTGTby/l-e^ 



^ jb{a)a'^ da 
y/r^ — a^t^ 



(2) 



where is the bulge M/L ratio, jbifl) the luminosity 
density at distance a from the center in the equatorial 
plane, and e = ^/l^^W/a? the constant intrinsic eccen- 
tricity. To determine jb from the brightness distribution, 
when e ^ we cannot invert the Abel equation (e.g., Bin- 
ney & Tremaine 1987] ) appropriate for spherical distribu- 
tions. An alternative approach (see Appendix) considers 
the "strip brightness", S'(a), defined as the integral of the 
bulge brightness distribution along a path orthogonal to 
the line of nodes, at distance a from the center. Under the 
hypothesis of transparency, this quantity is independent 
of the inclination of the galaxy. Since 



dS „ /- 
— = -27rVl 
da 



'ajb(a) 



(3) 



it is straightforward to determine jb (and hence the ro- 
tation speed from Eq. ^) once S{a) has been evaluated 
from the brightness distribution. Equation |^, in turn, can 
be expressed either in terms of a parametric or np pro- 
file respectively, depending on the kind of decomposition 
considered. 

2.3.2. The disk 

In the case of an infinitely thin exponential disk, the ro- 
tation speed is given by (Freeman 1970| ): 



vl{r)^4^Gn'^)rDV^[h{v)Ko{y)-h{y)Ki{y)] , (4) 
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Table 1. Galaxy sample and RC references. 



where y — r/{2rd), In and Kn arc the modified Bessel 
functions of first and second kind, respectively; S(0) — 
Tdl{0) is the central disk surface density, and the disk 
scale length. 

For the np disk we have adopted the following expres- 



sion for the circular speed (Toomre 1963; Kent 1986) 



vi{r) = 2nrG aY,{a) H{a,r) da 
Jo 

where S(a) is again the disk surface density, and 



(5) 



H{a,r) = 

irar 



K{k) - 

2E{k) 
7rr2(l-A:2)' 



k — -, for r < a (6) 
a 



for r > a 



(7) 



Here K{k) and E{k) are the complete elliptic integrals 
of first and second kind, respectively. Equation || can be 
rewritten in the more condensed form 



where pa is the only parameter of the distribution. This 
distribution, with Vc linearly increasing with r, is suited in 
some cases to represent the inner part of the halo RC; it is 
in fact the limit of Eq. ^ for r ^ r/^. We have adopted it 
whenever measurements of circular speed did not extend 
beyond the optical radius. In most cases the two models 
are mutually exclusive, since the latter cannot fit a flat RC, 
whereas Vh is not constrained if the data are restricted to 
the rising part of the RC. 

2.4- Modified Newtonian dynamics 

To test the predictions of modified Newtonian dynam- 
ics (M ONO: Milgrom |l983| ; Sanders |l990| ; Begeman et al. 



1988 - hereafter BBS), we have altered our model RC's 



according to the prescriptions given by BBS. In partic- 
ular, the relation between Newtonian (gn) and modified 
(g) acceleration 

gn gM(g/ao) 



v^{r) = AGr 



1 tE{t) 











Y.{rt)dt + 



(8) 



which we have used in our computations. 
2.3.3. The dark halo 

We have considered two different spherical distributions 
of matter to model the dark component. The first is a 
pseudo- isothermal sphere (Kent p.986 ) with density profile: 



with 

M(a;) =a;(l + x2)-i/2 

can be used to derive the modified expression for the cir- 
cular velocity. 

Although in most cases the available RC's are not 
very extended, in six cases (namely NGC 1024, NGC 2841, 
NGC 3593, NGC 4698, NGC 5879, and IC 724), the criti- 
cal acceleration (ao ~ 10^* cms^^) jg achieved within the 
radius sampled by the RC. These objects therefore can 
provide a test for MOND. 



p{r) 



47rG'(r2 



rl) 



and circular speed: 



(^) 



arctan 



(9) 



(10) 



where vt, the asymptotic circular speed, and rt, the scale 
length, are parameters of the distribution. The second is 
a sphere of constant density poi with 



4r2 



37rGpo 



(11) 



3. Results 

3.1. Rotation curve fitting 

With the exception of NGC 2841 and perhaps NGC 4698, 
our rotation curves do not extend far enough to constrain 
effectively both the stellar M/L's and the halo parameters 
of Eq. and result in high uncertainties on the param- 
eters (of order 50%) when all of them are fit simultane- 
ously. Moreover, a global fit of the three components, at 
least in the case of the isothermal halos, tends to attribute 
most of the mass to the halo at all radii. These difficulties 



were also noted by Kent (1988), who considered further 
assumptions and simplifications to better constrain the 
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parameter space. In particular he suggested three possi- 
ble kinds of fits. 

a) Maximum bulge-|-disk solutions (MBD hereafter), 
in which the maximum amount of mass compatible with 
the rotation observed is ascribed to the visible matter. 
This can be achieved by fitting the RC with only the two 
stellar components for radii within two disk scale lengths, 
or up to the disk velocity peak for the np case. When 
the observed RC is constant or rising beyond this limit 
the dark halo contribution is added. In our case the MBD 
approach can be justified by noting that these are all rel- 
atively high- luminosity objects with L ~ L*, so that DM 
should not be important in the inner regions (Kormendy 



1990| ; Salucci et al. |1991| ); in general the rising part of the 



RC is well reproduced, confirming that this hypothesis is 
probably correct. Moreover, in the case of NGC 2841 the 
simultaneous fit of the three components does not differ 
appreciably from the MBD one. For NGC 4698 it yields 
M/L ratios for disk and bulge lower by about a factor of 
two, although with higher uncertainties which make the 
result again consistent with the MBD results. 

b) Fits with fixed asympotic halo velocity, Vh- If we 
assume that this parameter can be determined from the 
flat portion of the RC, the number of free parameters is 
reduced to three. Only a few of our galaxies have extended 
RC's, so that an independent estimate of Vh is possible 
only for NGC 2841 and NGC 4698. In these two cases, 
however, the fitting routine is able to estimate Vh with 
reasonable accuracy (< 15%), so that we did not consider 
this kind of solutions. 

c) Fits with constant-density halos. Adopting the halo 
distribution described by Eq. |ri|, again we have only three 
free parameters. We find that in general such a halo does 
not affect appreciably the contribution of bulge and disk 
to the rising part of the RC, thus yielding M/L's consistent 
with the MBD values. 

On the basis of the above statements, and to define a 
homogeneous set of parameters, we decided to consider 
only the MBD solutions, assuming that the associated 
errors are realistic estimates of the true parameter un- 
certainties. The dark halo contribution was added for six 
galaxies out of 14; five of these have RC's measured at 21 
cm, while for the sixth, NGC 2639, the optical RC was 
sufficient to constrain the halo contribution. In the case of 
NGC 3593, the radio measurements, extending to about 
1.5 i?25, could be fit reasonably well without any dark 
component. For each galaxy we chose the model halo (the 
constant density or the pseudo-isothermal sphere) which 
yielded the best-quality RC fit; the constant-density halo 
was adopted for NGC 2639, NGC 4450, and NGC 5879. 
Comparing quantities which can be derived in both cases, 
such as the central density or the mass within the optical 
radius, we find no systematic difference between the two 
models. 



The resulting M/L's and halo parameters are given in 
Table |^, and Table || reports the masses of the various com- 
ponents. The values presented in the tables are the mean 
of the parametric and np results; the error is the largest 
of the values for the semi-difference (parametric and np) 
and the formal error on the fitted parameter. Figure ^ il- 
lustrates the fits to the observed rotation curves for both 
the parametric and np surface brightness decompositions. 
We note that the best results in most cases are obtained 
with the np fits, which in general match more closely the 
features of the observed curves. Quantitatively, comparing 
the obtained in the two cases, we find a median ratio 
of the parametric to the np value of 1.3. 

3.2. M/L ratios 

Disks have a mean M/L ratio of 1.6 and 1.0 (solar units) 
in J and K, with a dispersion of 0.6 and 0.4, respectively; 
the mean bulge M/L is 1.1 ± 0.6 in J and 0.6 ± 0.2 in K. 
We note that five galaxies (namely NGC 2639, NGC 3593, 
NGC 3898, NGC 4419 and NGC 5879) yield low values 
of bulge M/L; these were excluded from the mean calcu- 
lation (see Sect. 4.1). For NGC 4698, this happens only 



in the np case. Note that NGC 3593 and NGC 4419 are 
actively forming stars (Paper I), so a low M/L might 
be expected. Kent (1988) also found low M/L ratios for 
these five galaxies, and interpreted this as an inconsistency 
between photometric and kinematical data, possibly due 
to non-circular motions of the ionized gas in the inner 
part of the galaxies. The same conclusion was attained by 
Fillmore et al. ( 1986| ), who have modeled both rotation 
and velocity dispersions from emission and absorption-line 
data for six galaxies, three of which are also in our sample 
(see Table |l|). Their models suggest that in most cases even 
the observed emission-line velocity underestimates the ac- 
tual RC of the galaxy in the inner regions (r < 1 kpc). 
In the case of NGC 3898 and NGC 4450, for instance, a 
rough estimate of the bulge M/L from their model RC 
turns out to be respectively six and two times higher than 
ours. Even these 'corrected' values, however, remain lower 
by at least 15% than the disk M/L's estimated from the 
observed rotation. 



3.3. Dark halos 

The dark halo parameters we obtain are similar to the 
ones derived for other samples, in particular for later type 
spirals using photometry in the optical passbands (e.g., 
Kent 1986; Kent 1987). Scale lengths rh are comparable to 
the optical size of the galaxies and range from 18 to 30 kpc; 
the average halo mass is about 10" Mq within R25. The 
central halo densities show a rather narrow distribution, 



peaked at 30 -10 g cm with a dispersion of about 
20%. The average ratio of luminous to dark matter within 
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Table 3. Masses of the components. 



Fig. 1. Observed and fitted rotation curves. Upper and lower panels show respectively the results from parametric and np 
decompositions {K band). The dots are the observed curve; the dashed, dotted, and dot-dashed lines are respectively the 
contributions of bulge, disk, and dark halo to the model RC (continuous line). The mark on the abscissa corresponds to -R25. 



the optical radius is around 2, a typical value for bright 
spirals (PSS). 

3.4. MONO 



For the six objects selected in Sect. 2.4, we first performed 
a fit of the RC with qq and the bulge and disk M/L as 
free parameters. This yielded for ag a weighted mean of 
1.3 in units of 10~* cms~^. We then performed two sets 
of fits keeping qq fixed to our value of 1.3 and to the value 
of 0.8, corresponding to the ag = 1.21 obtained by BBS 
rescaled to Hq — 50 kms^^. 

A visual inspection of the observed RC's and the mod- 
els does not clearly favor one paradigm over the other. 
This qualitatitve statement is confirmed by the similar 
values of x^- Two representative cases are shown in Fig. 
|: NGC 2841, for which MONO yields a particularly good 
fit, and NGC 5879 for which a dark halo seems to produce 
a significantly better result. The choice of qq turns out to 
be crucial only for NGC 2841 where only ao = 1.3 yields 
a good fit. BBS also found that this galaxy required a 
higher critical acceleration, probably because of an error 
in the estimated distance, and did not include it in their 
average. After removing NGC 2841 from the sample, the 
remaining five galaxies yield ao — 0.74, fully consistent 
with the BBS value. 

4. Discussion 

Stellar content 

We have compared colors and M/L's for our sample with 
models of stellar population synthesis (SPS) to place con- 



straints on age and metallicity of the average stellar con- 
tent of bulges and disks. In particular we have consid- 
ered Worthey's model (1994, W94 hereafter) for single 
stellar populations with ages from 1.5 to 17 Gyr and dif- 
ferent metallicities, and the 1995 release of the model by 



Bruzual and Chariot (1993, BC95) for populations with 
solar metallicity and different star formation histories. 

The mean colors and M/L's of bulges and disks are 
plotted in Fig. || together with the two SPS models. The 
average M/L for the bulges are computed excluding val- 
ues below 0.2 in all bands. W94 colors and M/L's agree 
rather well with the observed values, whereas BC95 have 
both bluer colors and lower M/L's. The discrepancies be- 
tween the two models have been discussed by Chariot et 
al. (199(;). Inspection of the left panel of Fig. ^ shows 
that bulges are notably redder than disks, both in r — K 
and in J — K (see also Paper I). Nevertheless, the dis- 
crepancies between different SPS models, the degeneracy 
age/metallicity, and the possible effect of extinction make 
the color-color plot in Fig. ^ ambiguous as a diagnostic for 
distinguishing different stellar populations. 

Fortunately, M/L ratios disentangle, at least partly, 
the ambiguity of age and metallicity. We find that bulges, 
on average, have lower M/L's than disks in all bandail. 
Moreover, the values of M/L seen in the right panel of Fig. 
||, according to the SPS predictions, suggest that bulges 
are younger and more metal rich than disks. We note that 
both W94 and BC95 models of a given abundance follow 
approximately the same trend with age. Hence, the dis- 



^ J is not shown in the Fig. ^, but J-band M/L's behave in 
the same way as those in r and K. 
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Fig. 2. Fit with MOND to the RC's of NGC 2841 and NGC 5879 from parametric (upper panels) and np decompositions (lower 
panels). Dots represent the observed RC, whereas the continuous line is the model one. In both cases we plot the results for 
ao = 1.3. 



Fig. 3. Left panel: Comparison of average component colors to models of stellar population synthesis. The three thick lines with 
marked dots correspond to models (W94) with different metallicities as labelled: [Fe/H]=0.5 [Fe/H]=0.0, and [Fe/II]=-0.25. Age 
increases from bottom to top; designated points correspond to ages of 1.5, 2, 3, 5, 8, 12, 17 Gyr. The three thin continuous lines 
correspond to models (BC95) with fixed solar metallicity and different IMF's. Dots mark the same ages as for W94. models. 
Filled symbols represent the average colors of disks with the relative uncertainties; the triangle is for the parametric results and 
the square for the non parametric ones. Empty symbols are the bulge values. Right panel: the same comparison for M/L's (r 
and K bands). 



placement of the bulge and disk values relative to that 
trend implies that bulges are characterized by a younger 
age than that of the disks, independently of discrepancies 
between models. Except for extreme inclinations, i > 75°, 
dust in the central r egion s affects the bulge more than the 
disk (Bianchi et al. 1996 ). Consequently, the noted differ- 
ence between bulge and disk M/L cannot be attributed 
to internal extinction, since the correction would only in- 
crease the difference between the two. The possibility that 
the disk M/L's have been overestimated seems rather un- 
likely, even under the MBD hypothesis (although see Bot- 
tema|l997; Courteau & Rix 1997) since an average error 
greater than 30% on the estimated disk RC contribution 
would be required to make the disks' average age com- 
parable to bulges'. Finally, the bulge M/L's could have 
been systematically underestimated. This would be the 
case only if all the RC's (and not just the few already noted 
and excluded from the mean) rose too slowly in the inner 
regions with respect to the true circular velocity of the 
galaxy. We discussed briefly this point in Sect. 3.2, con- 
cluding that most likely such underestimates, if present, 
are not sufRcient to eliminate the observed difference in 
M/L between the components. 

That bulges may be more m etal rich than di sks is not a 
new resuh (Bica fc Al loin |l987| ; Dehsle & Hardy |l99l| ; Gio- 
vanardi & Hunt 199(\ , Paper I) , and abundance variations 
are thought to be driven by variations with mass (e.g., 
Zaritsky et al. 1994). That bulges appear to be younger 
than disks is somewhat surprising; the comparison with 
SPS models shown in Fig. ^ implies an age difference of 
around 50%, or 5 Gyr. Nevertheless, such a result may be 
interpreted in light of recent observational and theoretical 



work on bulge dynamics. Many bulges show kinematic and 
photometric signatures usually associated with disks, in- 
cluding flattened distributions, exponential fall-off, domi- 
nance of the rotation velocity component, and sp iral st ruc- 
ture in the bulge-dominated region (Kormendy 1993 and 
references therein). Moreover, some bulges have blue col- 
ors, the result of extremely young populations (Schweizer 



1990), and as noted in Paper I, at least three of the galax- 
ies in our sample appear to be actively forming starsB. As 
suggested by Kormendy and others, "bulges" may be built 
up over time from disk material transported to the central 
regions by gravitational perturbations; such bulges would 
appear younger than the disks from which they derive. 



4-2. Correlations with mass-to-light ratios 

Since, as for the photometric properties discussed in Pa- 
per I, disk characteristics are more reliably determined 
than those of bulges, we will concentrate on the M/L's 
obtained for the disks. 

Sever al authors (see for instance Djorgovski & Santi- 
ago 1993 , and references therein) have demonstrated that 
elliptical galaxies follow a relation which can be expressed 
in terms of a power law: M/L cx L" with a = 0.2 ~ 0.3 
depending on the sample and the photometric band. Kent 
(1986) also found that the disks of spirals follow a similar 
relation in the r band at fixed morphological type with an 
exponent a — 0.18 ± 0.07. More recently Persic, Salucci, 
and collaborators (PSS; Salucci et al. 1991) found for a 
sample of late-type spirals that in B, a ~ 0.30 ~ 0.35. 



^ Two of these have, however, been excluded from the calcu- 
lation of the mean. 
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Burstein et ah (1997) have suggested that a relation of this 
kind between M/L and luminosity is common to all the 
self-gravitating structures in the universe, ranging from 
globular clusters to clusters of galaxies. Finally, models 
of cold dark matter halos, based on N-body simulations 
and adiabatic infall for disk formation, require a variation 
of disk M/L with B luminosity in order to accommodate 
observed rotation curves (Navarro et al. 199(:). 

We have investigated the compatibility of the values of 
a found in the optical with our K-hand data, assuming the 
trend of M/L with L is due to disk stars and not to DM. 
We can convert the index as found in the B band to the 
K-hand value using the well-esta blishe d color -lum inosity 
relati on for spirals (Visvanathan 1981 ; Wyse 1982 ; TuUy 
et al. 1982 ). Based on recent data, Gavazzi ( 1993 ) finds: 



B -H 



-0.2 Mh + const. 



for early-type spirals. If we assume that H ~ K is inde- 
pendent of luminosity (which is likely since H — K is typ- 
ically small, ^ 0.2 mag), then we have a similar relation 
tor B- K with 

B - K = -0.2 Mr + const. . 

We take as =0.35 (PSS) together with the slope of the 
color-luminosity relation defined above, we infer a value of 
aK — 0.15. If we /it our data to a regression of M/L versus 
disk luminosity, we obtain aK — 0.22 ± 0.15, consistent 
with the expected value of 0.15. Figure ^ shows iC-band 
disk luminosity plotted against disk M/L (K), together 
with a line having slope a — 0.15. 

The M/L vs luminosity relation for elliptical galaxies is 
implicit in a more general relation, namely the one defin- 
ing the fundamental plane (FP) of ellipt ical galaxies (see 
review by Kormendy & Djorgovski |1989|): 



b 

c 



(12) 



where Re and /g are respectively the effective radius and 
surface brightness and CTc is the observed central veloc- 
ity dispersionEl. We have fitted a similar relation between 
disk scale lengths, central brightnesses, and peak rota- 
tion velocities of our disks, and find: a = 1.2 ± 0.2 and 
b — — 0.84± 0.1. This implies that the disks of early-type 
spirals also define a plane similar to that for ellipticals; 
the elliptical FP has a = 1.4 ± 0.15 and b = -0.9 ± 0.1. 
A similar result, but for the photometric properties only. 



We note that for elliptical galaxies the central observed ve- 
locity dispersion depends on total mass, whereas RC fitting in 
principle disentangles the contributi ons of stars and DM. If we 
assume, following Salucci & Persic 1997, that the fractional 



DM content is a decreasing function of luminosity also for el- 
lipticals, an M/L vs L relation for the stellar component alone 
would require a steeper slope than the one derived from the 
FP relation. 



was reported in Paper I, both for bulges and disks. When 
we derive the M/L vs L relation from Eq. |l^ and the virial 
theorem we obtain: 



M/L cx L^I{0) — 

which contains a residual dependence on the central 
brightness I{0). However our values for a and b yield 

M/L cx L"-3/(0)-°-°7 

where, as for ellipticals, the dominant dependence is on 
luminosity. 

As discussed by Djorgovski & Santiago (1993), a re- 



lation between M/L and luminosity (or mass) can come 
about in several ways. One possibility is that the disk 
M/L's are contaminated by a DM contribution which has 
a density profile similar to that of the stellar disk. The 
sense of the M/L vs L relation would require this DM 
fraction to increase with luminosity, contrary to the trend 
observed for the global DM fraction which increases with 
decreasing luminosity (e.g., PSS). Alternatively, the MBD 
hypothesis could be incorrect, and the stellar disk M/L 
constant. In this case, though, the trend in Fig. ^ would 
require the MBD hypothesis to be more valid in lower 
luminosity systems, contrary to common beliefs. 

Another possibility is that disks of different luminosi- 
ties harbor different stellar populations, which is also sug- 
gested by the color-luminosity relation mentioned above. 
W94 predicts that at fixed age, initial mass function 
(IMF), and star formation rate (SFR), M/L is an increas- 
ing function of metallicity in the optical, but a decreasing 
one in the NIR. This suggests that the M/L vs L correla- 
tion cannot be understood in terms of a metallicity varia- 
tion. Alternatively, such a correlation could be driven by a 
change of average age or star formation history with lumi- 
nosity. Again, the observed difference in the slope of the 
correlation at different wavelengths can be compared to 
the predictions of SPS models. To test this possibility we 
made use of the BC 95 m odels, at fixed (solar) metallicity 
and IMF (Salpeter |l955| ) , considering single burst popu- 
lations at different age T, and populations with different 
e-folding time t (exponential SFR) at fixed age (10 Gyr). 
We have approximated the model dependence of M/L on 
T and t with power laws, whose index in B and K has 
been determined with a best fit. Assuming M/L^ cx L^-^^ , 
in the first case the model predicts a difference in the 
slope of this relation of 0.17 passing to the K band, con- 
sistent with the predicted difference of ~ 0.2. Considering 
the variation with r the prediction is ~ 0.2, in similarly 
good agreement. It, therefore, seems plausible that the 
variation of M/L with luminosity is driven by different 
star formation histories, and the consequent different stel- 
lar mixes. 
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Fig. 4. The M/L vs luminosity for the disks {K band) with the trend derived from as and the color- luminosity relation 
described in the text. 



4-3. Dark halos 

We claim to recognize the presence of a dark halo in six 
of our galaxies: three of them are pseudo-isothermal and 
three constant-density spheres. We do not find any sys- 
tematic difference between the two models, at least in 
terms of central densities or masses within R25 . The statis- 
tics are sparse, but we can attempt a comparison with the 
general relations found by PSS in the B band for a large 
sample of late- type spirals. In particular, we can check 
that the ratios of the halo central density to the critical 
density, po/pc, and the dark to visible mass at the optical 
radius, Mh/My, are consistent with the correlations with 
the B-band galaxy luminosity given in PSS. They find: 



Mv 



= 0.4 



Lb 

T * 



-0.9 



and ^ = 3.5.104f^ 

Pc 



-0.07 



with L*, 



101"-4l|. 



If we assume the same color- 



luminosity relation as in Sect. 4.2, we find in the K band 
after scaling to iJo = 50 kms^^: 



2) NIR M/L's appear to vary with disk luminosity, 
consistently with the trend in the B band, and the B — K 
color-luminosity relation. Such a trend can be probably 
imputed to a systematic variation of the disk average stel- 
lar population with luminosity, due to a variation of either 
star formation history or age, but not of metallicity. 

3) We find no clear evidence to prefer either standard 
Newtonian gravitation plus dark halos or modified gravi- 
tation, as the quality of the fits under both paradigms is 
comparable. The value of the critical acceleration param- 
eter oo that best accommodates the sample as a whole is 
1.3 • 10~** cms-2. 

4) The dark halo parameters we obtain for six of our 
galaxies are roughly in agreement with the values obtained 
for late type spirals from _B-band photometry. In partic- 
ular our data are consistent with the correlation, already 
noted by PSS, between dark-to-visible mass ratio and lu- 
minosity. 



Mv 



-0.7 



= 0.8 — and 



Lk 
L 



K 



^=3.7. 

Pc 



10^ 



T * 



-0.06 



with Pc = 3H§/8t:G = 0.5 • lO^^g g ^m-^ and = 
—23.97, the absolute magnitude corresponding to L*j^. We 
estimated M}^ from the color-luminosity relation, with 
the zero order coefficient fixed by the median values of 
A4b and B — K for our sample. 

The trends shown in Fig. |^ are consistent with the an- 
ticorrelation between dark and luminous mass found in 
PSS. Moreover, there is no striking discrepancy between 
our galaxies and the behavior of later type systems, sug- 
gesting that dark halos are similar for all spirals. In the 
right panel of Fig. ^ our data reveal roughly the same Mh 
as found in late-type galaxies with the same B luminosity 
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5. Summary 

We summarize here the main conclusions of this paper. 

1) The comparison of bulge and disk NIR average col- 
ors and mass-to- light ratios to stellar population synthesis 
models suggests that the bulges of our sample galaxies are 
younger and more metal rich than the disks. 
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Fig. 5. Central halo density and ratio of dark to visible mass vs galaxy K absolute magnitude. Filled triangles correspond to 
isothermal halos, open circles to constant density ones. The dashed lines represent the NIR relations between these parameters 
corresponding to the ones obtained by PSS in the B band. 



A. The strip brightness method to derive the lu- 
minosity density of an axisymmetric elHpsoid 

The strip density method was first introduced by 
Schwarzschild (1954) to evaluate the gravitational poten- 
tial of the Coma cluster. In the following we describe the 
method and extend it to the more general situation of 
an oblate, axisymmetric, ellipsoidal distribution of mass 
(and luminosity) . This technique offers a quite simple and 
straightforward way to recover the spatial density distri- 
bution, and the potential and rotation curve of such sys- 
tems, especially when they are not described parametri- 
cally but rather observed as projected distributions on the 
sky 

Let us consider an ellipsoidal distribution, centered at 
O, with oblate rotational symmetry about z, inclined by 
an angle i to the line of sight [Fig. ^(b)]. The isodensity 
surfaces can be parametrized by their semimajor axis, that 
is their equatorial radius a: 



j{x,y,z) = j{a) , 
with a 



1 -e2 



(Al) 



where the constant e = 1 ^ ^ is the intrinsic eccen- 
tricity. The projected distribution on the plane of the sky 
(^, 77) has elliptical symmetry with an apparent eccentric- 
ity e sin i. The surface brightness is then 



y 



2 t 2 



l-e2 

and area A(a,Co) — tta/ 1 — (a^ — ^o^) 



(A4) 



Then, 



j{x = Co,y,z) dydz 



j{a)dA = 2ti\J\ - e2 / j{a)ada (A5) 

"'Co 



and 



(A6) 



In practice, the observed image is divided in strips nor- 
mal to C and integrated in the 77 direction. This gives <5'(Co) 
for Co spanning from to the outermost radius of the im- 
age. Equation is then used to derive j(a), and with it 
the potential and circular speed. From a numerical point 
of view the method is easy to implement and was tested 
on a number of analytical solutions. Due to the strip in- 
tegration prior to the differentiation in Eq. A6, the noise 
introduced is low and hardly appreciable even in the outer 
regions. 



n^,v) = Hi) , 



with q = W ^2 



1 — (e sin z)2 



(A2) 



We now define the strip brightness at distance Co from 
the image minor axis [Fig. ^(a)] as 



^(Co) = 



-|-oo 



(A3) 



On the other hand 5 (Co) is the luminosity of a section 
perpendicular to a; at x = Co of the ellipsoid. Such a sec- 
tion will have again elliptical symmetry with eccentricity 
e. Any elliptical isodensity contour of this section can be 
described by: 
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Fig. 6. Geometry of the ellipsoid (a) as seen projected on the sky plane, and (b) from a lateral viewpoint. The symmetry axis 
z is inclined with respect to the line of sight and the axes x and ^ coincide with the line of nodes, interception of the equatorial 
plane {x^y) with the sky plane {£,,rj). The strip brightness is the integral of the surface brightness along a strip perpendicular 
to the line of nodes, as sketched by the dashed lines in panel (a). 
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Table 1. Sample galaxies and RC references. 



Name 


Distance 


R25 


Bt 


RC 






[Mpc] 


[arcsec] 


[mag] 




(1) 




(2) 


(3) 


(4) 


(5) 


NGC 


1024 


69.2 


117 


13.08 


1 


NGC 


2639 


70.4 


55 


12.56 


1 


NGC 


2775 


27.1 


128 


11.03 


1 


NGC 


2841 


19.2 


244 


10.09 


2 


NGC 


3593 


17.2 


157 


11.86 


1,3,4,5 


NGC 


3898 


34.8 


131 


11.60 


1,6,7 


NGC 


4378 


51.4 


87 


12.63 


1 


NGC 


4419 


23.9 


99 


12.08 


1 


NGC 


4450 


23.9 


157 


10.90 


6,8 


NGC 


4698 


zo.y 


110 


11 Afi 


1,5, y 


NGC 


4845 


27.1 


150 


12.10 


1 


NGC 


5879 


26.2 


125 


12.22 


6,10 


NGC 


6314 


139.2 


43 


13.80 


1 


IC 


724 


125.6 


70 


13.40 


1 



Notes to Table 1. 

Column 2: see paper I. Columns 3-4'- values from the 
Third Reference Catalog of Bright Galaxies (de Vau- 
couleurs et al. 1991) Column 5: References for the 
RC's (O = optical, R = radio). 1: Rubin et al. 1985 
(O): 2: Begeman 1987 (R); 3: Bertola et al. 1996 (O); 
4: Demoulin 1969 (O): 5: Krumm & Salpeter 1979 
(R); 6: Fillmore, Boroson & Dressier 1986 (O). For 
NGC 4450 we considered just the rising part of their 
observed curve, since the declining outer part is incon- 
sistent with the 21 cm data from reference 8; 7: van 
Driel & van Woerden 1991 (R); 8: Guhathakurta et 
al. 1988 (R); 9: Warmels 1986 (R); 10: Broeils & van 
Woerden 1994 (R). 



Table 2. Fits to the rotation curves: mass-to-light ratios and halo parameters. 



Name n {M/Lj)i {M/Lj)d (M/L^),, (M/L^)d po Vh 

[IQ-^^g an-^] [kpc] [km s^^] 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 


1024 


3 










0.81 


(0.20) 


1.20 


(0.45) 








NGC 


2639 


4 


0.00 


(0.44) 


1.58 


(0.25) 


0.06 


(0.06) 


0.91 


(0.15) 


43.3 (9.4) 






NGC 


2775 


3 


1.39 


(0.20) 


1.51 


(0.15) 


0.80 


(0.10) 


0.86 


(0.08) 








NGC 


2841 


2 


1.44 


(1.27) 


2.31 


(0.27) 


0.76 


(0.70) 


1.21 


(0.14) 


26.4 (5.5) 


21.4 (2.7) 


309 (10) 


NGC 


3593 


1 


0.15 


(0.05) 


0.74 


(0.15) 


0.06 


(0.03) 


0.42 


(0.14) 








NGC 


3898 


2 


0.14 


(0.10) 


1.96 


(0.11) 


0.06 


(0.06) 


1.22 


(0.15) 


42.4 (10.0) 


18.7 (7.0) 


368 (50) 


NGC 


4378 


3 


0.78 


(0.30) 


2.51 


(0.48) 


0.52 


(0.21) 


1.55 


(0.31) 








NGC 


4419 


3 


0.16 


(0.09) 


1.27 


(0.23) 


0.07 


(0.05) 


0.72 


(0.20) 








NGC 


4450 


3 


0.50 


(0.09) 


1.49 


(0.37) 


0.34 


(0.01) 


0.81 


(0.18) 


28.4 (7.7) 






NGC 


4698 


2 


0.32 


(0.31) 


2.63 


(0.57) 


0.19 


(0.19) 


1.97 


(0.77) 


31.6 (19.5) 


30.7 (13.0) 


335 (50) 


NGC 


4845 


4 


1.88 


(0.33) 


1.39 


(0.18) 


0.73 


(0.14) 


0.77 


(0.11) 








NGC 


5879 


3 


0.08 


(0.08) 


0.89 


(0.30) 


0.05 


(0.05) 


0.58 


(0.19) 


19.1 (0.8) 






NGC 


6314 


2 


0.85 


(0.13) 


1.31 


(0.19) 


0.38 


(0.05) 


0.82 


(0.15) 








IC 


724 


2 


1.57 


(0.30) 


1.71 


(0.19) 


0.70 


(0.11) 


0.98 


(0.06) 









Columns 2-5: M/L's are in solar units 



Notes to Table 2. 



Table 3. Masses of the components. 



Name 


Mt 


Md 


Mv 




Mtot 


Mh/Mv 


(1 


) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


NGC 


1024 


17.26 (0.47) 


9.74 (1.35) 


27.01 




27.01 




NGC 


2639 


0.86 (0.40) 


25.24 (2.74) 


25.46 


18.54 (3.59) 


44.00 


0.73 


NGC 


2775 


8.62 (1.63) 


18.35 (2.76) 


26.97 




26.97 




NGC 


2841 


5.91 (0.53) 


28.65 (4.54) 


34.57 


11.67 (1.52) 


46.24 


0.34 


NGC 


3593 


0.12 (0.08) 


1.71 (0.22) 


1.83 




1.83 




NGC 


3898 


0.81 (0.84) 


13.23 (2.02) 


14.04 


17.48 (1.30) 


31.52 


1.24 


NGC 


4378 


5.29 (1.06) 


22.60 (1.88) 


27.89 




27.89 




NGC 


4419 


0.24 (0.14) 


4.91 (0.72) 


5.16 




5.16 




NGC 


4450 


1.25 (0.49) 


13.05 (2.64) 


14.30 


10.60 (2.90) 


24.90 


0.75 


NGC 


4698 


1.69 (0.13) 


15.74 (4.75) 


16.59 


3.96 (2.10) 


20.55 


0.24 


NGC 


4845 


2.67 (4.55) 


7.51 (2.06) 


10.17 




10.17 




NGC 


5879 


0.02 (0.02) 


1.69 (0.15) 


1.71 


4.74 (0.19) 


6.45 


2.77 


NGC 


6314 


9.81 (1.39) 


20.83 (2.07) 


30.64 




30.64 




IC 


724 


20.19 (4.61) 


30.05 (3.79) 


50.24 




50.24 





Notes to Table 3. 

Column 2: bulge mass; Column 3: disk mass; Column 4- bulge + disk mass; Column 
5: halo mass; Column 6: total mass (all values are in units of 10^° Mq). Column 7: 
dark-to-visible mass ratio. 



1,0 - 



u 



o 
to 



0.0 - 



-1.0 - 



1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



J I I I I I I I I I I I I I I I I L 



-23 -24 -25 

Absolute luminosity (K-mag) 




-23 -34 -25 -36 
Absolute luminosity (ff-mag) 



in 



J I I I I I I I I I I I I I I I I I I I L 



0- % 



_ v. 



0.0 h 



-1.0 - 



-S.O 



578 



1 



I I I I I I I I I I IV I 



-33 -34 -35 -36 
Absolute Imnlnoslty (ff-mag) 



